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ABSTRACT 

It  is  shown  that  the  space  of  bivariate  C1  piecewise  cubic  functions  on 
a  hexagonal  mesh  of  size  h  approximates  to  certain  smooth  functions  only  to 
within  0(h3)  even  though  it  contains  a  local  partition  of  every  cubic 
polynomial. 
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SIGNIFICANCE  AND  EXPLANATION 


One  measures  the  approximation  power  of  a  family  S  of  piecewise 
polynomial  approximating  functions  on  some  partition  in  terms  of  the 
meshsize  h  of  that  partition.  Typically,  the  error  of  approximation  goes  to 
zero  like  hr  as  the  meshsize  goes  to  zero,  with  r  depending  on  the 
smoothness  of  the  function  being  approximated.  There  is  a  maximal  r  typical 
for  the  space  s  used,  and  faster  convergence  rates  are  possible  only  for 
very  special  functions.  Naturally,  one  would  like  this  optimal  rate  or 
approximation  order  hr  to  be  as  fast  as  possible,  i.e.,  would  like  the 
maximal  r  to  be  as  large  as  possible.  In  any  case,  it  is  an  important 
practical  question  to  ascertain,  for  a  given  approximating  space  S,  what  its 
optimal  approximation  order  is.  ^ 

In  a  multivariate  context,  this  is  not  an  easy  question  to  answer,  even 
for  pp  spaces  on  very  simple  partitions,  as  soon  as  one  demands  that  the 
approximating  functions  have  a  certain  amount  of  smoothness.  The  report  shows 
that  a  certain  line  of  attack  which  was  thought  to  be  quite  reasonable  by 
several  workers  in  the  field  (including  the  authors)  cannot  succeed,  in 
general.  This  makes  the  problem  even  harder  than  it  was  initially  thought  to 
be. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
suaswury  lies  with  KRC,  and  not  with  the  authors  of  this  report. 


Approximation  order  from  bivariate  C1 -cubical  A  countaraxampla 
C.  M  Soor  and  X.  Hflllig+ 


•  We  dial  with  a  soali  (d^)  of  approximating  spaces,  generated  fro*  a 


fixad  apaoi  I  by  a  a lip la  aoallngi 


■h  *"  Vf)  * 


(ohf)(x)  i-  f(x/h)  ,  all  f,  x,  h  . 


Wa  ara  lntaraatad  In  tha 


obtalnabla  fro*  ,  i.e.,  In 


dist(f,  Sn) 


aa  a  function  of  h  and  for  f  sufficiently  smooth.  Hers, 

dlattf,  Bjj)  i-  Inf  If  -  gl 
and  1*1  la  tha  sup  nor*  on  ao*a  oloaod  doaaln  o  tr  , 


If  I  i“  aup _ Iftacl  I  . 


It  Is  easy  to  see  that 


dlattf ,8b)  -  o(*£(h>)  ,  all  f  e  C(g) 


In  oas*  8  contains  a  local  and  stable  partition  of  unity.  By  this  we  Man  that 


■  on  8? 


for  som  ^  e  S  with 


aup^  die*  supp  ♦,  <  • 


,upl  '♦i*  * 


The  last  condition  is  autoaaticelly  satisfied  in  case  the  4^  are  all  nonnagatlva. 
W*  are  interested  in  suitable  conditions  on  8  which  insure  that 

dlsttf,  8b )  ■  0(hk)  for  all  sufficiently  smooth  f  . 

It  is  easy  to  see  that  (1)  Implies 


This  aaterial  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  Mo.  MCS-7927062,  Mod.  1. 
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*k  £■  *  ' 

with 

«•  polynomials  of  degree  <  k  . 

On  the  other  hand,  thia  oondition  is  clearly  not  sufficient  for  (1)  since,  e.g. ,  8  - 

implies  that  Sj,  ■  for  all  h  ,  hence  dist(f,  8^)  is  independent  of  h  in  this  case. 
What  seeas  to  be  needed  is  that  P^  be  contained  in  8  "locally",  such  as  is  con¬ 
tained  locally  in  8  in  case  8  contains  a  local  partition  of  unity. 

Here  is  a  precise  formulation  of  such  a  condition. 


Condition  **  .  For  every  p  €  8^  ,  there  exists  a  ’ sequence'  <f. )  in  8  so  that 
P  "  and  supi  dian  supp  ^  . 


At  least  in  case  8  is  "uniforn",  e.g.,  8  is  a  space  of  piecewise  polynomial 
functions  on  some  uniforn  partition  of  H2  ,  one  would  expect  to  conclude  from  this  the 
validity  of  (1),  i.e.,  global  approximation  order  0(hk)  (see,  e.g.,  [BD],  ID]). 

It  is  the  purpose  of  thia  note  to  give  an  exaepla  of  a  piecewise  polynomial  space  8 
on  a  uniform  partition  of  H2  which  satisfies  Oondition  P4  ,  yet  gives 

dlst(f,  8h)  >  constf  h3  (2) 

for  some  positive  constf  and  for  the  particular  function 

f«x»— 9  (x<1)x(2))2  . 


this  dashes  all  hopes  that  the  approximation  order  from  a  piecewise  polynomial  scale 
(8h)  could  be  settled  by  finding  out  which  polynomials  are  contained  locally  in 
8  .  Presumably,  scam  stability  has  to  be  added  to  Oondition  before  global  approximation 
order  0(hk)  can  be  deduced. 

Here  is  an  outline  of  this  note.  In  section  2,  we  show  that  cur  example  space  8 
satisfies  Oondition  P^  and  that  (8^)  has  approximation  order  0(h3)  at  least.  In  Section 
3,  we  identify  8  within  a  larger  space  of  piecewise  cubic  functions  as  the  annihilator  of 
a  set  A  of  local  linear  functionals.  We  also  show  that  there  exists  a  bounded  ad  ** \ 0 
so  that  a(X)  If  -  0  for  all  f  with  oonptct  support,  while  I _ v_„ 


-2- 


a(X)A  »  0 


only  If  tho  sum  la,  In  effect,  ovar  all  or  over  non*  of  A  .  Thla  la  important  for  tha 
final  a notion,  in  which  wo  prova  that  (B^)  haa  approximation  order  at  moat  0(h3)  . 

2.  anooth  pp  fumotioma  and  bom  apl Ine*.  No  oonaidar  bivariate  pp  (<■  plaoewlee  poly¬ 
nomial)  function*  on  tha  partition  i  of  ^  obtained  from  tha  thraa  famlllaa  of  maah- 


x(1)  -  n  ,  x(2)  an,  x(l)  -  x<2>  +  n  ,  all  n  «  S  . 

we  ara  particularly  interested  in  tha  apaoa 


»  -  *4, A  -  Va"®1^ 


of  piecewise  cubic  functiona  on  tha  partition  A  and  in  C1  . 

Ho  have  forogon*  tha  opportunity  to  make  tha  aymnatrlaa  in  A  mora  apparant  by  hmkg 

o 

tho  thraa  famlllaa  of  maahlinaa  intaraact  aaeh  othar  at  an  angle  of  120  (aa  la  dona,  e.g. , 
in  [Fr] ) .  Thla  naadlaaaly  complioataa  tha  notation.  It  la  aufficlant  to  not*  that  any 
permutation  of  tha  maahlina  famlllaa  can  be  accomplished  by  aoma  linear  map  on  B?  ,  and 
tha  corraapondlng  change  of  variables  leaves  P*  4  invariant. 

\  a table  and  local  partition  of  unity  In  8  la  constructed  In  [80],  [BHjJ  aa  follow*. 
Conaider  tha  linear  nap  Pi B?  — >  a?  character laad  by  the  fact  that 

pi  • 

P*i  -  <  *2  '  i-2'5  /  ' 

Ul4**  '  1-5  J 

with  the  ith  unit  vector  (in  **) .  Let  N  be  tha  P-ahadow  of 

B  »-  [0, 1) 5  , 

i.a.,  M  la  the  diatribution  given  by  tha  rule 

M*  I-  /  AOP  ,  all  p  . 

B 

Since  M  la  the  ahadow  of  a  box,  we  call  it  a  bom  aplina.  It  la  immediate  from  tha  defin¬ 
ition  that  N  >  0  ,  aupp  N  -  P(B)  ,  and 

S  !•(•—))  *  1  ,  with  J  t-  I2  , 

J« 

the  laat  bacauaa 

«(-!)♦  "  *°P  “  ll  3  •"»-/.♦• 

3  3  B+3  r*[0,1J3  ** 


Further,  on*  verifies  that  MSS. 

In  addition,  [BHj]  providas  an  l^ -bounded  linaar  functional  X  with  support  in 
aupp  N  ao  that 

p  «  E  Xp(*+j)  M<— j>  ,  for  all  p  8  P  (1) 

jej 

and  show*  how  this  result  leads,  in  standard  quasi-lnterpolant  fashion,  to  the  conclusion 
that 

disttf,  Sj^  -  0(h3)  (2) 

for  all  sufficiently  saooth  f  .  On  the  other  hand,  [BHj]  makes  clear  that  (1)  is  sharp. 


i.e.,  that 

»4  \  8M  #  0  , 

with 

8,,  «-  span  (n<  —  3>)jeJ  • 

Therefore,  (2)  provides  the  optimal  approximation  order  from  (SM>h)  . 

Since  S||  is  a  proper  subspace  of  8  ,  this  only  provides  a  lower  bound  on  the 
approximation  order  from  (8h)  .  Further,  8  satisfies  Condition  »4  .  To  see  this,  recall 
from  [M2]  that  S„  contains  certain  cubic  polynomials,  e.g.,  the  two  polynomials 

xt—>  x(2)2(3x(1)  -  x(2) )  and  xw  x(1)J(3x(2)  -  x(1))  .  (3) 

In  addition,  the  linear  change  of  variables  which  Interchanges  a,  and  e,+e2  ,  hence 
oarries  e2  to  -e2  ,  leaves  A  invariant,  hence  leaves  8  invariant.  Therefore,  8 
contains  the  cwoic  polynomials  which  result  from  such  a  change  of  variables  from  those  in 
(3).  Up  to  parabolic  terms,  these  are  the  polynomials 

x»— »  (x(  1)-x(2) )2(2x(  1)  ♦  x(2))  and  x»— 0  x(1)2(2x(1)  -  3x(2)),  (4) 

and  the  four  polynomiala  (3),  (4)  are  linearly  independent  since  the  matrix 


"0  0  3  -1  " 

-13  0  0 

2-301 
.  2  -3  0  0  . 


of  their  leading  ooeffieiants  ia  invertible 


J#  Imttia  ooogdi— t— >  in  this  auxiliary  auction,  wu  identify  B  ■  » '  as  a 

4, A 

aubepaee  of  P®  A  satisfying  carta In  homogeneous  conditions.  We  find  it  aoat  convenient  to 
express  thaaa  oonditiona  in  tana  of  tha  Barnatain  ooordinataa  for  pp  funotiona  on  a 
triangulation,  aa  introduced  by  farin  [F],  following  aarliar  work  by  do  Caataljau  tC]  and 
Sabin  (■] .  Sara  ia  a  abort  explanation  of  thin  vary  uaafui  roproaantation. 

on  a  aingla  triangle  T  with  wartiooa  0,  V,  and  W  ,  wo  uae  barycantria  ooordinataa. 
Thia  aaana  that  each  point  x  ia  aaaooiatad  with  tha  tripla  («,t,w)  for  which 

x  ■  uU  +  W  ww  ,  and  u  +  y  +  w  "  1  . 

In  thaaa  tana,  wa  describe  a  polynomial  p  of  degree  (n  by 


i+j+k-n  liX  *i3X  * 


4  (X)  ■  - — -  uiylw,t  . 

fijk'  '  it jlkl  u  w 

Wa  deal  with  tha  31  choices  in  thia  representation  by  associating  b^*  with  tha  point 

xijk  **  (iO  *  tf  *  *W>/n 
for  all  i  ♦  J  +  k  -  n  .  Tha  resulting  function 


bt  xljk  bijk 

ia  called  tha  a(arnatain  or  -tsisr )-net  for  p  (with  respect  to  t  ).  It  ia  Independent  of 
how  wa  associate  tha  vertices  of  t  with  tha  letters  B,  V,  w  .  Moreover,  if  a  is  an 
affine  change  of  variables,  than  tha  B-net  for  pOA  (with  respect  to  a'^t)  )  is  bon  . 

This  aakes  it  easy  to  compare  polynomial  pieces  across  triangle  adgas.  For  example,  on 
an  edge  of  t  ,  p  is  entirely  determined  by  b  raatrietad  to  that  edge.  Moreover,  if 


P’  is  a  polynomial  of  degree  <  n  on  a  triangle  t’  having  that  edge  in 


with  T  , 


then  p  ■  p*  on  that  edge  iff  b  «  b*  on  that  edge,  with  b*  the  »-net  fox  p*  . 

Higher  smoothness  across  such  an  edge  is  also  very  simply  expressible  in  terms  of  b 
and  b'  ( see  [f] ) .  We  now  describe  these  conditions  only  to  the  extent  that  we  need  them, 
i.e.,  we  describe  the  conditions  which  an  f  «  p®  ^  must  satisfy  to  belong  to  s  *  pj  ^  . 
Since  such  an  f  ia  continuous,  tho  B-nets  of  its  various  pisoes  must  agree  at  all  points 
of  overlap  in  their  domains.  We  can  therefore  think  of  the  B-nets  of  its  various  pieces  as 
forming  one  B-net,  a  function  b{  defined  on  all  of 


*&"!?> 


J3  i-  (M/3)2  .  (1 

Let  t  ,  t*  b«  two  triangles  of  A  with  *  eo— on  edge  e  .  Thor*  or*  four  points  of  J3 
on  c  .  Bach  of  tha  thraa  paira  Xj,  x2  of  aueh  noighboring  points  has  a  nsarsst  Jj- 
naighbor  y  in  T  and  a  nsarsst  Jj-naighbor  y 1  in  T'  and  togathsr  thsss  foor  points 
font  tha  vertices  of  a  parallelogram  halvsd  by  c  .  Ons  nay  verify  dirsotly  that  f  has 
continuous  first  dsrivativss  across  c  if  and  only  if 

bjtx,)  +  bf(x2)  ■  bf(y)  ♦  bf(y*) 
for  sach  of  thsss  thraa  para  Holograms. 

Thus,  assooiatsd  with  saoh  sdgs  c  in  A  thsrs  ara  thraa  linaar  functionals 
.  0 

A  on  w,  .  ,  of  tha  font 

If  i-  bg(Xj )  ♦  b^Xj)  -  bf(y)  -  bf(y')  , 
with  Xj  noighboring  Jj-points  on  c  and  y  ,  y*  adjacant  Jj-polnta  in  tha 

noighboring  triangles.  Mots  that  wa  have  so  normalised  X  that  tha  edge  points  receive 
weight  1  and  tha  off-adga  points  receive  weight  -1  . 

Since  A  contains  three  distinct  edge  types,  this  gives  altogether  nine 
nonoverlapping  alasees  ,  i,  j-1 ,2,3,  of  linear  functionals.  To  be  precise,  we 

associate  A^  with  segMnt  i  of  edge  j  ,  and,  in  particular,  A^  with  the  middle 
segment.  He  need  not  be  more  precise  then  that. 


Figure  i.  The  nine  classes  of  local  linear  functionals 


each  class  A^  la  la  ft  Invariant  undar  tranalatlon  of  tha  Independent  varlabla  by 


aoaw  j  a  J  •  Thalr  union 


A  A„  U...U  AJ3 


charactariaaa  S  within  r.  .  in  tha  aanaa  that 

A*4 

a  “  ker  A  i-  n  k«r  X  . 

.  XflA 

Next,  wa  aaak  t  f  r\0  ao  that 

Z  a(X)  X  ■  0  . 

xeA 

By  thia  wo  naan  that 

I  a(X)  Xf  -  o 

for  ovary  f  of  compact  aupport.  For  such  an  f  ,  tha  aua  haa  only  finitaly  many  nonaaro 

tana  ainca  aaoh  X  haa  aupport  only  in  aoaa  pair  of  adjaeant  triangles  of  A  . 

Thera  may  bo  many  solutions,  but,  when  wa  require  additionally  that  a  bo  constant  on 

each  A^  ,  than  tha  solution  sat  can  be  shown  to  be  four-dimonaional.  There  is  a  two- 

dimensional  aat  of  solutions  whioh  vanish  on  ,  all  j  .  These  solutions  are  of  no 

interest  to  us  since  for  them  wa  have  already 

Z  a (X)  X  -  0  (2) 

suppXcG 

for  various  bounded  sets  6  . 

Solutions  a  for  which  (2)  only  holds  if  the  sum  ia,  in  effect,  over  all  or  over  none 
of  A  are  obtained  as  follows. 


Let 


at/ 


be  such  that  a 


I A 


A  - 


A(i,J)  , 

with 

-a 

-8 

-Y 

2a 

28 

2Y 

—a 

-8 

-Y 

and  a  +  8  Y 


Then  Z  a(X)  X  ■  o  . 


Proof.  By  its  definition,  each  X  carries  an  f  to  a  weighted  aua  of  values  of  bf 
We  can,  therefore,  understand  Ea(X)X  by  ooaputing  the  weight  it  assigns  to  bf(x)  for 
each  x  «  J3  .  There  are  three  types  of  points  x  in  J3  ,  those  at  a  vertex  of  A  , 


•7' 


those  inside  an  edge,  and  those  inside  a  triangle.  Ns  consider  each  type  in  turn. 


A  vertex  of  A  serves  as  an  edge  point  for  six  X's  ,  one  from  each  of  the  classes 
A^  with  i  ?  2  .  Thus,  the  total  weight  at  a  vertex  point  is 

Ei*2Mi,j)  -  -2 ( o+B+Y )  -  0  . 

As  to  inside  edge  points,  consider  without  loss  of  generality  one  on  an  edge  of  type 
1.  Such  a  point  serves  as  an  edge  point  for  one  1  e  A21  and  one  X  e  u  ,  and  as 
an  off-edge  point  for  one  X  e  A^  uAJ2  and  one  X  e  u  AJ3  .  Its  total  weight  is 
therefore 

A(2,1)  +  A( 1 , 1 )  -  A( 1 ,2 )  -  All, 3)  - 

2a  +  (-a)  -  (-8)  -  (-Y)  -  o  ♦  8  +  Y  *  0. 

Finally,  an  interior  point  is  an  off-edge  point  for  three  X's  ,  one  each  from  A^  , 
j-1,2,3.  Its  weight  is  therefore 

-20  -  28  -  2Y  -  0  .||| 


We  extend  each  X  to  the  continuous  linear  functional  Xi  on  Cl*2)  with  the  aid  of 
the  local  linear  map  I  which  associates  f  with  the  unique  element  If  of  P°  A  which 
agrees  with  f  on  J3  .  It  is  then  a  simple  matter  to  check  that,  for  every  f  6  *5  ,  the 

map  X  •— >  Xif  is  constant  on  each  A^  (since  each  X  vanishes  on  *4  and  ,  for  any  j 

f  and  f(*+j)  differ  only  by  some  cubic  polynomial).  In  particular,  for 

fix)  i-  Ixl 1 ) x ( 2 ) ) 2  , 

we  have  t using  the  association  of  edga  types  Indicated  in  Figure  1)  that 


(«|A  ) 
'  ij 


Therefore, 


-6-6  6 

3  3  -12 

-6-6  6  J 


£  a(AAj)  Uf)lA4il)  -  18<o  8)  -  36Y 


i,j-1 


U 


The  number  *  is  nonsero  for  many  choices  of  o,  8,  Y  for  which  o  +  8  +  Y  "  0  .  Hake 

such  a  choice.  Then,  for  a  square  Q  with  sides  parallel  to  the  axes, 

£  »(X)  Xif  ■  x*area(Q)  +  01 perimeter (Q ) ) 

•uppXcg 

while  *,Upp*cg  *>*s  eupport  only  on  triangles  near  the  boundary  of  g  ,  hence 


(3) 


(4) 

I  £  a (X)  \1  I  •  0( perimeter (Q) )  . 

suppicQ 

This  is  the  essential  fact  required  in  the  next  and  final  section. 


4.  An  nppnr  bound  for  the  approximation  order.  In  this  section,  we  establish  the  main 

point  of  this  note.  With  S  ■  p]  .  as  defined  earlier,  we  show  that,  for  all  small 

4,C 

enough  h  , 

dist(f,  S^)  >  const  h3 

for  some  positive  const  and  for  the  particular  function 

fixt— ■>  (x(1)x(2))2  . 

For  the  proof,  we  pick  some  axis-oriented  square  Q  in  G  and  consider 

wh  1  supp\5yh  a(D  I  °1/,h 

with 


ah  i“  {  x/h  i  x  e  Q  }  . 


Since  Sh  £_  ker  li^  ,  we  have 


Further, 


dist(f,  Sh)  >  distQ(f,  Sh)  >  distfi(f,  ker  uh>  *  lvhf l/,Wh,c(Q) 


Vh£  -  h  f  ' 


while  area(Q^)  ■  area(Q)/h2  and  perimeter(Qh)  ■  perimeter (Q)/h  •  Therefore,  from  (3.3), 

lMj.fl  "  const  h2  +  Oth3) 
n 

with  const  i“  | <  I  area(Q)  >  0  ,  while,  from  (3.4), 


This  shows  that 


•Veto)  ■  0(1/h) 


disttf,  Sjj)  >  const  h3  +  Oth4) 


for  some  positive  const  ,  as  asserted. 
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